Conformal thin-sandwich puncture initial data for boosted black holes 
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We apply the puncture approach to conformal thin-sandwich black-hole initial data. We solve 
numerically the conformal thin-sandwich puncture (CTSP) equations for a single black hole with 
non-zero linear momentum. We show that conformally flat solutions for a boosted black hole have 
the same maximum gravitational radiation content as the corresponding Bowen-York solution in the 
conformal transverse-traceless decomposition. We find that the physical properties of these data are 
independent of the free slicing parameter. 
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I. INTRODUCTION 

Numerical simulations of binary black-hole collisions 
require initial data. A black-hole collision scenario of 
current interest involves two black holes in orbit before 
they collide. As the black holes orbit, they emit gravita- 
tional radiation, and their orbits' radii slowly decrease. 
Eventually the two black holes are so close that no sta- 
ble orbits exist, and they plunge together and merge. 
Gravitational radiation emission circularizes the orbits, 
so we expect the final orbits before merger to be almost 
circular. For many reasons long-term numerical evolu- 
tions are difficult, and therefore ideal initial data for a 
binary black-hole collision simulation would describe two 
black holes in an almost circular, slowly decaying orbit 
just prior to their merger, as well as all of the gravita- 
tional radiation that has built up during their slow in- 
spiral. We do not yet know how to construct such initial 
data. However, a number of techniques exist to construct 
initial data for two boosted black holes, and to identify 
initial-data sets that exhibit some of the characteristics 
of circular orbits 0, 0, IE 0, 0, El ■ I n this paper we 
outline a new procedure to construct initial data for mul- 
tiple boosted black holes, and present solutions for single 
boosted black-hole spacetimes. 

A popular starting point for constructing binary black- 
hole initial data is the Bowen-York solution of the mo- 
mentum constraint 0, E|. This solution can be used 
to construct two black holes with arbitrary prescribed 
masses, linear momenta and spins. An effective-potential 
technique was developed by Cook 0] to choose values 
of these parameters that correspond to quasi-circular or- 
bits. This method also allowed a prediction of the loca- 
tion of the innermost stable circular orbit (1SCO) of the 
two black holes. 

In Cook's original approach, the black holes were con- 
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structed on a two-sheeted topology, where the sheets were 
connected at the black- hole throats. The initial data was 
inversion symmetric between the two sheets. This ap- 
proach allowed inner boundary conditions to be placed 
on the black-hole throats, and the remaining initial-value 
equation (the Hamiltonian constraint) was solved on only 
one sheet, with the region inside the throats (the second 
sheet) excised from the computational domain, avoiding 
all coordinate singularities. A later approach by Brandt 
and Brugmann removed analytically the coordinate 
singularity at the location (or "puncture" ) of each black 
hole, and made it possible to solve the Hamiltonian con- 
straint without excising any region from the computa- 
tional domain, and without any need for inner bound- 
ary conditions. Baumgarte [l^] used Cook's effective- 
potential method to construct initial data for binary 
black holes in quasi-circular orbits using the puncture 
approach. The predictions of the ISCO for inversion- 
symmetric and puncture Bowen-York data were in good 
agreement. 

The Bowen-York solution arises from making certain 
choices of the freely specifiable quantities in the con- 
formal transverse-traceless (CTT) decomposition of Ein- 
stein's equations These choices are made because 
they make the equations easier to solve, and not for any 
strong physical reasons. As such, we have cause to doubt 
that these data are astrophysically realistic and we seek 
alternatives to Bowen-York data that are more physically 
motivated. 

A recent alternative decomposition of the initial-value 
equations is York's conformal thin-sandwich (CTS) de- 
composition 01 . The free quantities in the CTS de- 
composition are more closely linked to dynamics than 
those in the CTT decomposition. In particular, we can 
make choices consistent with a quasi-equilibrium space- 
time, which is what we expect for two black holes in 
quasi-circular orbit . 

Binary black-hole initial-data sets have been con- 
structed using the CTS decomposition by a number 
of groups. Grandclement, Gourgoulhon and Bonazzola 
0, 01 solved a variant of the CTS equations using exci- 
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sion techniques, to construct conformally flat initial data. 
Cook and Pfeiffer 0, 13 a l so use d excision techniques, 
but with improved inner boundary conditions. Yo, et. 
al. used a simplified version of the inner boundary 
conditions suggested in |l5j. 

It would be useful to be able to solve the CTS equa- 
tions in the puncture approach 0, Without the 
need for complex inner boundary conditions, the punc- 
ture approach is potentially easier to implement than 
excision. CTS-puncture binary black-hole initial data 
would provide a useful comparison with excision data 
results, and also allow evolution codes that employ punc- 
tures ISISIlOl to take advantage of CTS data. 

The puncture approach was extended to the CTS equa- 
tions in 01 ■ There it was pointed out that it will be dif- 
ficult (and perhaps impossible) to construct CTS binary 
black-holes in quasiequilibrium in the puncture frame- 
work. It will be important to investigate the impact of 
this within quasi-circular orbits in CTS-puncture binary 
data, but we will not address that issue in this paper. For 
the moment we will confine ourselves to outlining a gen- 
eral procedure for constructing multiple boosted black- 
hole spacetimes, and apply our method to single boosted 
black holes. We will leave the question of orbits for a 
later publication. 

Our starting point is a result by Laguna L 2lj, who 
found analytic solutions to an approximation to the con- 
formal thin-sandwich puncture (CTSP) momentum con- 
straint. These solutions describe a single conformally-flat 
black hole with non-zero linear or angular momentum, 
and reproduce the Bowen-York form of the conformal ex- 
trinsic curvature. Laguna's single-hole solutions provide 
an analytic result with which to test our numerical code, 
and also motivate the technique we employ for the full 
coupled CTSP system. This technique, which amounts 
to specifying a value for the shift vector at the punc- 
tures (s), allows us to construct initial-data sets for one 
or more black holes, each with non-zero linear momen- 
tum. 

In Section ^ we review the conformal thin-sandwich 
decomposition, and in section II I II outline its incorpora- 
tion into the puncture framework. We consider Laguna's 
solution for a single black hole with non-zero linear mo- 
mentum in Section llVl We describe our numerical meth- 
ods in Section We present numerical solutions of the 
full CTSP system for a single boosted black hole in Sec- 
tional an d consider their physical properties. 



timeslice. In the 3+1 decomposition the spacetime met- 
ric can be written 



ds 2 



-N 2 dt 2 + jij (dx l + (3 l dt) {dx 1 + (i ] dt). (1) 



The lapse function N gives the proper time between 
slices, Ndt, and the shift vector f3 % describes the coor- 
dinate drift between slices. Using these quantities one 
can calculate the extrinsic curvature using 
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{dtjij - Vifo - V,f3 t ) 



(2) 



The constraint equations on one slice are, in vacuum, 
R + K 2 - KijK tj = 0, 



0. 



(3) 
(4) 



The covariant derivative V, is taken with respect to the 
spatial metric, as is the Ricci scalar R. The quantity K 
is the trace of the extrinsic curvature. 

The initial data, 7y and Kij , evolve in time according 
to the evolution equations 



d t K i:j 



-2NK l3 + Vipj + Vjpi, (5) 
-ViVjiV + N (Rij - 2K lk K) + KK t j) 
+l3 k V k K t , + K ik VjP k + K kj ViP k . (6) 



We use the conformal thin-sandwich decomposition of 
the initial- value equations [TsT ]. The spatial metric 7y is 
related to a background (or conformal) spatial metric 7y 



7ii 



(7) 



The extrinsic curvature is split into its trace and trace- 
free parts, 



Ki 



A, 



1 



K. 



(8) 



The trace of the extrinsic curvature K is the same in 
both the physical and conformal space, K = K. We 
provide a conformal weighting for the lapse function, N — 
i/j 6 N, and introduce the time derivative of the conformal 
metric, 



dt7ij, 



(9) 



II. CONFORMAL THIN-SANDWICH 
DECOMPOSITION 



In the ADM (3+1) decomposition of Einstein's equa- 
tions [H, |22, the ten Einstein equations are projected 
onto spacelike hypersurfaces, resulting in four constraint 
equations that the metric, 7^, and extrinsic curvature, 
Kij, of each timeslice must satisfy, and six evolution 



which transforms as Uij = ij) 1 uy. 

In this decomposition the conformal extrinsic curva- 
ture is given by 



X-- — 

2N 



(10) 



and its conformal weighting Ay = ip 2 ''Aij and A 1 -? 



equations that tell us how 7^ and evolve to the next -0 W A V follows. The resulting conformal thin-sandwich 
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decomposition of the initial-value equations is, 
uum, 



1 



5r/2 



~2* K 



V, 



2N 



' 4 ■ A % 3 



1 - 4 



J V27V 



(11) 
(12) 



V*(ip''N) = Nil)' ( li; 



' 4 ■ 4 lJ 



+■0" 



where 



A L /3* 



= v*^' + - -f J Vfc/3 fc , 
3 J 



(13) 

(14) 
(15) 



Equation {TTJ is the CTS version of the Hamiltonian con- 
straint @ ■ The momentum constraint |@J is now written 
as l(12|). Equation 113|l follows from specifying <9 t /T in the 
ADM evolution equation for 7\T, obtained from the trace 
of ijHJl. This equation provides us with a lapse function, 
AT = tp 6 N, and is necessary to complete the CTS system 

M 

In this decomposition the free data are the conformal 
metric 7-y and trace of the extrinsic curvature K, and 
their time derivatives on the initial slice, and d t K. 
Throughout this paper we make the choices of conformal 
flatness, 7^ = (the flat metric in Cartesian coordi- 
nates), maximal slicing, K = 0, and initial stationarity 
of K and the conformal spatial metric, Uij — 8 t K = 0. 
These stationarity choices are useful in the context of 
binary black holes in quasi-circular orbits 0, E3 • With 
these choices, the conformal thin-sandwich equations be- 
come 

1 



A L /3* - OL/S^'VjbiJV = 0, 

v 2 (v> 7 ao = AV' 7 



7^ 



' 4 ■ - 4*-? 



(16) 
(17) 
(18) 



Equation l|18|l is the maximal-slicing equation. It is 
a condition that, if the current slice is maximal, and we 
evolve to the next slice with a lapse that obeys (T8J), then 
that slice will also be maximal. 



With 



0, the conformal extrinsic curvature 110(1 



is now given by 



Ai 



2N 



(19) 



III. THE PUNCTURE METHOD 



We wish to solve the equations 1(16( 1 - ((18(1 using the 
puncture approach [Tl| . In this approach the conformal 



factor for n black holes is written as 



x - rrij 



(20) 



Here m, are parameters that characterize the mass 
of each black hole. In the special case of a single 
Schwarzschild black hole, mi equals the ADM mass of 
the black hole. The rj are the distances of each black hole 
from the origin. Equation ((21) (I gives the Brill-Lindquist 
conformal factor for n momentarily stationary black holes 
[281, 129I] . plus a function u, which describes the deviation 
in the conformal factor due to a non-zero extrinsic curva- 
ture. The function u is regular over all of the conformal 
space [ll| . The topology of this solution consists of one 
hypersurface that contains n black holes, and each black 
hole connects through an Einstein-Rosen bridge to an ad- 
ditional hypersurface, giving a total of n+l hypersurfaces 
connected by n Einstein-Rosen bridges. 

In the CTSP approach the conformal lapse N is 

also split into an analytic singular and unknown regular 
part, 



Nip 7 = NiP = 1 + J2tt + v 



(21) 



We are free to choose the constants c,; they determine 
the value of the lapse at spatial infinity on the other hy- 
persurfaces. However, negative choices of Ci will lead to a 
lapse function that goes through zero, and this will cause 
division-by-zero problems in a numerical construction of 
Aij via (|T§1) . For this reason, we choose positive values 
of Ci in this paper. We investigate the effect of different 
choices of a on the physical results in section IvTl 

With these choices, the conformal thin-sandwich punc- 
ture (CTSP) system is 



V 2 u 



V 2 v 



1 



A L /3 l -(L/3)^V,lnA> 



Nil) 7 



0. 



-J-*A %3 A^ 



(22) 

,(23) 
(24) 



Brandt and Brugmann 11] showed that the solution u 
of ((22|) will be C 2 everywhere on R 3 if the conformal 
extrinsic curvature A^ diverges no faster that 1 /r f at 
each puncture. The same will be true for 1(23|) IT]]- In 
the construction of Aij via ((19(1 . the conformal lapse N 
goes to zero as rf. If Aij is to satisfy the divergence 
requirements of ((22(1 and l(23|) . then (L/3)y must go to 
zero at least as fast as rf. If this requirement is met, 
then the CTSP equations will be regular over all of R 3 , 
eliminating the need for any region of the computational 
domain to be excised. 

The particular solution that we obtain will depend on 
the boundary conditions we apply, both at the outer 
boundary and at the puncture (s). There is no need for 
puncture boundary conditions on equations 1(22(1 and 1(231) 
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if their source terms are smooth at the punctures. On 
the other hand, puncture conditions are needed for the 
momentum constraint l|24|l. if we wish to avoid a trivial 
solution. For example, consider a single boosted black 
hole. We expect the shift vector to fall off to zero at spa- 
tial infinity. However, if we impose the outer boundary 
condition [3 l — on the shift vector, and do not place any 
condition on the shift vector at the puncture, a numerical 
code will yield the trivial solution j3 l = everywhere. A 
zero shift vector will lead to a zero extrinsic curvature, 
and will not represent a boosted black hole. 

In this work we will show that a black hole with non- 
zero linear momentum can be generated by imposing a 
value on the shift vector at the puncture, in addition to 
outer boundary conditions. The full CTSP system is 
solved at every point on a Cartesian grid, except at the 
punctures, where the momentum constraint is not solved, 
and instead at the ath puncture we require that 



# = 4- 



(25) 



The values A l a parametrize the linear momentum of the 
black hole, which will point in the direction of the vector 
A l a = A l a /\A\ for each hole. In this prescription the shift 
vector will be zero at spatial infinity. We will provide a 
motivation for this procedure in the next section, where 
we consider solutions for a single black hole with linear 
momentum. 



IV. A SINGLE BLACK HOLE WITH LINEAR 
MOMENTUM 

Laguna |2l| has found analytic perturbative solutions 
of the CTSP system ~ <E2) that describe a single 
boosted or spinning black hole. Laguna's solutions yield 
the Bowen-York extrinsic curvature. In the case of a 
boosted black hole, these solutions are perturbative to 
second order in the momentum for u and v, and linear 
for (3 l . We will focus on the momentum constraint, which 
takes the form Laguna considered if we choose u = v = 
(or u = v = 1 in Laguna's notation). The conformal 
factor and lapse splittings (|2U[) and l|21|l are then 



m 
2r 



NiP 7 = 1 + 



2r 



(26) 
(27) 



The lapse parameter c corresponds to the parameter —b 
in [2l| . The solution of the CTSP momentum constraint 
for a Bowen-York black hole with linear momentum P % — 
(P, 0, 0), located at the origin in Cartesian co-ordinates, 
is 



-Ax 2 f(r, to, c) — g(r, to, c) 
20(m + 2r) 6 

xyf(r,m,c) 
~5(TO + 2r) 6 ' 

xzf(r,m,c) 
' 5(m + 2r) 6 ' 



P 



(28) 



(29) 



(30) 



with 

f(r,m,c) — (m + c)m 2 + 12(to + c)mr 

+60(m + c)r 2 + 160r 3 , (31) 

g(r,m,c) = 5(5m + c)to 4 + 60(5m + c)m 3 r 
+2(749m+ U9c)m 2 r 2 
+8(497m + 97c)mr 3 
+120(49m + 9c)r 4 + 4480r 5 . (32) 

A calculation of the conformal extrinsic curvature via 
(PI yields 

A ij = i [PV + PV + (fi - n l n 3 )P k n k ] , (33) 

where n % = x % jr are normal vectors directed away from 
the puncture. Equation (|33|l is the Bowen-York extrinsic 
curvature for a black hole with linear momentum P % . 

As r — - > oo, the conformal lapse approaches unity, 
N — ► 1, and the momentum constraint l|24|) approaches 
its conformal transverse-traceless (CTT) decomposition 
form, 



A h W l = 0. 



(34) 



The solution W l that corresponds to a boosted Bowen- 
York black hole with momentum P l is 



W l = -— {IP 1 + n i n,P j ) 

The asymptotic form of - J3DJ is 

I3 l = 2W\ 



(35) 



(36) 



The solution - ESJ is not unique: ft + V i is also 
a solution of the momentum constraint, for any constant 
vector V. This choice will determine the value of the 
shift vector both at the puncture and asymptotically. 
The value of the shift vector at the puncture, /3q, is 



05 



5m - 



4to 2 

and its asymptotic form is 



-P l + V\ 



2W l + V 1 



(37) 



(38) 



In the solution presented in l|28|l - l|30|l . we have cho- 
sen V 1 = 0, and the shift vector is (3 l = [—(5m + 
c)P/(4m 2 ), 0, 0] at the puncture, and asymptotically ap- 
proaches twice the Bowen-York vector potential for a 
boosted black hole. Laguna's form of the solution is 
zero at the puncture, and has the asymptotic form (3 l — 
(5m + c)P i /(4:m 2 ) + 2W\ 

This solution motivates the condition on the shift vec- 
tor at the puncture that we introduced in section ITTT1 To 
numerically construct a single black hole with linear mo- 
mentum, one can solve (|22H - (|24p. requiring that, at the 
puncture, 



(5m 



4m 2 



-P l 



(39) 
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where c > 0, and apply an outer boundary condition con- 
sistent with the Bowen-York vector potential (|35[l . We 
will show in section [""/] that this technique does indeed 
produce a black hole with linear momentum, although 
the magnitude of the linear momentum will be equal to 
P only when we make the assumption that u = v = in 
the solution of the momentum constraint. Generalization 
to two (or more) black holes is straightforward. 

Given the CTSP system of equations (2J - we 
must choose appropriate outer boundary conditions. Re- 
quiring that our solutions be asymptotically flat implies 
that we can apply 1 jr Robin boundary conditions to the 
functions u and v. We are therefore assuming that be- 
yond the outer boundary of the numerical grid the func- 
tions u and v behave as, 



r 
fe 2 



(40) 
(41) 



where k\ and /c2 are constants. Note that in a numerical 
code the finite outer boundary will introduce errors, and 
as a result the values of k\ and &2 may vary across the 
surface of the outer boundary. To impose 1 [r fall-off on 
a function /, we impose 



N%(rf) = 0, 



(42) 



where N l is the unit normal to the outer boundary. The 
outer boundary condition l|42|) is applied to the solu- 
tions of both the Hamiltonian constraint l|22l) and the 
maximal-slicing equation (|23[l . 

In this paper we will use three types of outer boundary 
condition on the shift vector. In Section El we will use a 
Dirichlet outer boundary condition and prescribe a value 
for the shift vector at each point on the outer boundary, 



If = B l 



(43) 



where B l is a known vector. In the case where we choose 
v = 0, B l is the analytic Laguna solution 128|) - (|30[l . 



In the general CTSP case we could choose B l = 0, but 
this would introduce an outer boundary error that falls 
off only linearly as the outer boundary is moved out. 

The second type of outer boundary condition we use is 
a scalar Robin boundary condition, as described above 
for the solutions of the Hamiltonian constraint and 
maximal-slicing equation, and given by (|42|l . The er- 
ror in this boundary condition will fall off quadratically 
as the outer boundary is moved out. 

A third type of outer boundary condition makes use 
of our prior knowledge of the asymptotic angular de- 
pendence of the shift vector. As we pointed out earlier, 
the CTS momentum constraint approaches the form of 
the CTT momentum constraint for large r. The asymp- 
totic behavior of any solution of the CTT momentum 
constraint with symmetries consistent with a boosted 
black hole [3(| suggests that the shift vector will have 



the asymptotic form 



k. 



IP 1 + P k n k ri 



(44) 



where P l is a unit vector in the direction of the linear 
momentum. In the case of a black hole boosted in the re- 
direction, we choose P % = (1,0,0). The resulting Robin 
outer boundary condition is 



k 



N j d 3 (r(3 l ) = — [(N i - 27i l N ] 7ij)P k n k + n\P N j )] . 

Note that this is an approximate vector Robin bound- 
ary condition since we are assuming a direction for the 
momentum P l in using (|44|l to determine k% separately 
for each component of the shift vector. 

All three of the shift outer boundary conditions H42J) . 
(14*31 and (|4*5*|) are investigated in Sections M and IVTl 



V. NUMERICAL TECHNIQUES AND 
CONVERGENCE TESTS 

We solve the CTSP system - <G2l with a multi- 
grid elliptic solver, which is a modification of the 
BAM_Elliptic solver in the Cactus infrastructure |3lj . 
The CTSP equations are coded with a second-order 
finite-differencing stencil as a coupled elliptic system on 
a Cartesian grid. 

In all solutions presented here the code imposes a punc- 
ture condition of the form (|39|l . This is done by placing 
the black-hole puncture at a grid point on the numerical 
grid. All of the CTSP variables are well-behaved at the 
puncture, so there is no problem placing the puncture on 
a grid point. We require (|39[) at the puncture, and do 
not solve the momentum constraint there. The remain- 
ing two CTSP equations, (|22|l and l|23() . are solved at all 
points on the numerical grid. 

As a first illustration of these techniques, we consider 
the Laguna solution, i.e., we choose u = v = and solve 
only the momentum constraint. We employ Dirichlet 
outer boundary conditions, imposing values on the shift 
vector at the outer boundary equal to the solution 128|) 
- (|30[) . making the choice m = 1, c = 1 in the decompo- 
sitions of ip and Nip 7 in (|20|) and (|21|l . At the puncture 
we impose the condition (|39|l . We solve the momentum 
constraint for a single black hole with linear momentum 
in the x-direction of magnitude P — 0.5m. Convergence 
plots for this case are shown in Figures [*"] and [21 which 
represent solutions with outer boundaries of 2.0m and 
0.02m respectively. These outer boundaries would be far 
too close for a general problem with Robin outer bound- 
ary conditions on the shift vector, but are possible in 
this case, where the analytic outer boundary values are 
known. All solutions were found with grid sizes of 32 3 , 
64 3 and 128 3 points. The solution with outer boundary 
at 2m therefore has a base resolution of h = 0.125m. 

The figures show the disagreement between f3 l and the 
solution (|28|l - (|30|l for three grid resolutions. If the 
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solutions were second-order convergent, the errors would 
fall by a factor of four each time the grid resolution h was 
halved. First-order convergence would cause the errors 
to fall by a factor of two. 

The solutions in Figure^show worse than second-order 
convergence, but better than first-order convergence. We 
found that the convergence behavior slowly improved as 
finer resolutions were used, but even at the extremely fine 
resolution of ft. = 1.25 x 10~ 3 rn (the base resolution in 
Figure we do not see full second-order convergence. 
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FIG. 1: Convergence behavior of the errors in the momentum 
constraint solution along the z-axis (f3 x ) and the 3D diagonal 
{f3 v and /3 Z ), with outer boundary at 2m. 



The poor convergence behavior appears to be due to 
the use of the puncture condition i|39|) . This condition 



FIG. 2: Convergence behavior of the momentum constraint 
along the z-axis (f3 x ) and the 3D diagonal ((3 y and f5 z ), with 
outer boundary at 0.02m. 



effectively applies an inner boundary condition at a single 
point, unlike a true inner boundary condition, in which 
a condition is applied to an entire inner surface. 

The effect of the puncture can be seen by solving the 
system with a true inner boundary. A cubical region 
around the puncture was excised, and the analytic solu- 
tion imposed on the resulting inner boundary. Figure |3| 
shows the convergence behavior of the ir-component of 
the shift vector for this case, with an outer boundary at 
1.0m and a base resolution of h = 0.0625m. It is clear 
that in this case we have second-order convergence. The 
use of a puncture condition reduces the rate of conver- 
gence. This is unfortunate, but not catastrophic: we still 
have a convergent code, and since we know the correct 
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x=0y=0z=0 r — > do 

tt even even even 1/r 

v even even even 1/r 
(3 X even even even 1/r or (7 + i 2 /r 2 )/r 
fP odd odd even 1/r or xy/r 3 
f3 z odd even odd 1/r or xz/r 3 

TABLE I: Coordinate plane symmetries and outer boundary 
conditions for CTSP variables for a single boosted black hole 
with momentum in the ^-direction. 



condition to apply at the black-hole puncture, but do not 
know what inner boundary condition to apply on an ar- 
bitrary inner boundary, the use of a puncture condition 
remains advantageous. 
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FIG. 3: Convergence of [3 X (along z-axis) when an inner 
boundary condition is used, and the outer boundary is at 
1.0m. Differences were scaled assuming second-order conver- 
gence. 



VI. A SINGLE BOOSTED CTSP BLACK HOLE 

We are now ready to consider the full CTSP system, 
where the Hamiltonian constraint (|22[) and maximal- 
slicing equation l|23|) are coupled to the momentum con- 
straint (|24fl . In order to save computer resources, we 
solve the full CTSP system on one octant of the numer- 
ical grid (i.e., for only positive x, y and z) by employing 
the known symmetries of the solutions. In the case of a 
single black hole boosted in the x-direction, the CTSP 
variables exhibit the symmetries listed in Table [I] Table 
[I] also shows the Robin outer boundary conditions that 
we use. We will apply both scalar and vector Robin outer 
boundary conditions to the shift vector, as described in 
Section HVI and compare their effect on the solutions. 

Figure 0] shows a convergence plot for a single black 
hole with the outer boundary at 16m, and a base grid 
spacing of ft = 0.125m. For these solutions a scalar 
Robin outer boundary condition (|42|l was used for all 



CTSP variables, including the shift vector. The differ- 
ences between solutions with resolutions ft, 2ft/3, ft/2 
and ft/3 are plotted, corresponding to grid sizes of 128 3 , 
192 3 , 256 3 and 384 3 points. The differences are scaled 
assuming second-order convergence. These results are 
consistent with the test problem considered in section IVI 
the solutions are between first- and second-order conver- 
gent. In particular, (3 y and /3 Z show poor convergence 
near the puncture, and we can also see that convergence 
deteriorates a little near the outer boundary. 

Figure [5] shows convergence results for the same prob- 
lem, but with the vector Robin outer boundary condition 
(|45|l applied to the shift vector. These plots show that 
both types of Robin outer boundary conditions have the 
same convergence properties. 

Figure [5] shows the CTSP solutions for this problem, 
with ft = m/24. The solutions of u, v and (3 X are almost 
independent of whether we use scalar or vector Robin 
outer boundary conditions on the shift vector. However, 
the solutions of (3 y and (3 Z differ significantly at the outer 
boundary, and the solutions for both choices of shift vec- 
tor outer boundary condition are shown in Figure El 

Having generated solutions to the CTSP equations, we 
now calculate a number of physical quantities on the re- 
sulting spacetime slice. 

In the analytic solution l|28l) - (|3U|) we know the re- 
lationship between the value of the shift vector at the 
puncture and the linear momentum of the black hole. 
This relationship does not hold for the full CTSP sys- 
tem. We calculate the linear momentum of the initial 
data using 0, 0] 

P l = A^SSj. (46) 

The surface integral is computed at the outer boundary 
of the computational grid. The integral was constructed 
using global Killing vectors of the conformal space (which 
are asymptotic Killing vectors of the physical space) and 
can in fact be computed at any radius, so long as the 
surface surrounds the puncture. For the setup we have 
described, the symmetries will enforce P v = P z = 0. 

As stated above, the linear momentum of the solution 
is known analytically in the Laguna case, with ill'L'l) and 
(|23|1 solved using the Bowen-York extrinsic curvature in 
the source terms. This test case allows us to measure 
the accuracy of the solutions with respect to a physical 
quantity. Figure shows the error in the linear momen- 
tum for a single boosted black hole with m = 1, c = 1, 
Ph = (—0.75,0,0), and with the vector Robin boundary 
condition (|45|l applied to the shift vector. The analytic 
solution has P x = 0.5m. The code was run with two dif- 
ferent resolutions ft = 0.125m and ft = 0.0625m, and at 
three different choices of outer boundary location, 16m, 
24m and 32m. We find that the location of the outer 
boundary has a far greater effect on the accuracy of the 
solutions than the resolution. As such, for subsequent 
runs we choose ft = 0.125m and an outer boundary at 
32m as being an optimal combination of accuracy and 
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FIG. 4: Convergence of the full CTSP system for one boosted black hole with outer boundary at 16.0m, using scalar Robin 
outer boundary conditions on the shift vector. The differences are scaled assuming second-order convergence. The base grid 
spacing is h = 0.125m. 



economy of computational resources. 

Figurc[5]shows the linear momentum in the x-direction 
as a function of the puncture value of (3q for full CTSP 
solutions with vector Robin outer boundary conditions 
on the shift vector. Figure [5] shows how this value differs 
from that of the Laguna solution; as we would expect, 
the Laguna solution is a good approximation for small 
values of 

For a given initial-data set, we can also calculate the 
total ADM mass, Madm- In the conformal space, this is 



given by 

M ADM =-— I VVd 2 ^. (47) 

In terms of the puncture splitting of the conformal factor 
H20JI the ADM mass can be rewritten as a volume integral 

m 

n 

Madm = + — / ^A^A^dV, (48) 

i n J 

The ADM mass integral is over all space, but the nu- 
merical grid covers only a finite region. In order to esti- 
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FIG. 5: Convergence of the full CTSP system for one boosted black hole with outer boundary at 16m, using vector Robin outer 
boundary conditions on the shift vector. The differences are scaled assuming second-order convergence. The base grid spacing 
is ft = 0.125m. 



mate the contribution to the integral (|48|) from beyond 
the numerical grid, we need to estimate the values of the 
CTSP variables in that region. Following Baumgarte 01 
we do this by assuming that beyond the numerical grid 
the CTSP variables behave exactly as we required in their 
respective outer boundary conditions. For example, we 
required that 



(49) 



at the outer boundary. Given some radial vector from 
the origin to a point on the outer boundary of the grid, 
we can calculate the constant k\ at that point. The value 
of u at any point along that line, beyond the grid, can 



now be estimated. The accuracy of fci depends on the 
location of the outer boundary and the resolution of the 
grid, but the calculation of u beyond the grid is consistent 
with the outer boundary condition. 

A similar technique can be applied to the shift vector 
and v. However, an alternative approach is to recall that 
at the outer boundary the shift vector will behave like 
the Bowen-York vector potential for a single black hole, 
as described in the discussion before 144|) . When the con- 
formal extrinsic curvature is calculated using (|19f) it will 
resemble the Bowen-York extrinsic curvature (|33|l at the 
outer boundary, up to a multiplicative factor. That fac- 
tor will be related to the linear momentum of the black 
hole. Instead of calculating that factor numerically, we 
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FIG. 6: Solutions of the full CTSP system for one boosted black hole with outer boundary at 16m, with h — m/24. The 
solutions for u, v and (5 X did not differ substantially when either scalar or vector Robin outer boundary conditions were applied 
to the shift vector. The solutions of /3 y and /3 Z are shown for both types of outer boundary condition. 



can simply use the linear momentum calculated with 14(j|) 
and construct the extrinsic curvature outside the numer- 
ical grid with (|33|) and use this in the ADM-mass integral 
(|48|l . The values of (3 % and v outside the numerical grid 
are now unnecessary. The code we used to calculate the 
ADM mass integral was adapted from a code written by 
Miranda Dettwyler to perform the same calculation on 
Bowen-York puncture data |33| • 



Given the ADM mass and total linear momentum of a 
data set, we can estimate the maximum possible amount 
of gravitational radiation that these initial-data sets con- 
tain. The maximum gravitational radiation content can 



be defined as 

~W - M ' ( j 

where M is the puncture ADM mass of the black hole, 
defined as the ADM mass as calculated on the second 
hypersurface ^l| • The puncture ADM mass is also often 
referred to as the bare mass of the black hole 0, . 
The puncture mass for a single black hole is given by 

M = m{l + u ), (51) 

where uq is the value of the function u at the punc- 
ture and m is the mass parameter in (|20f) . Figure ITUI 
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FIG. 7: The percentage error in P x for different choices of 
resolution and outer boundary for the Laguna case, with c = 
m = 1, /ft = -0.75, and analytic P 1 = 0.5m. 



full CTSP system, but ignoring u and v in the construc- 
tion of the conformal extrinsic curvature. This is sim- 
ilar to the Laguna case described in Section |VJ except 
that now the Hamiltonian constraint and maximal-slicing 
equations are also solved, and u and v, although they are 
ignored in constructing the terms in the momentum con- 
straint, are not zero. All solutions were found with an 
outer boundary of 32m, with resolution h = 0.125m. Due 
to small systematic errors in the calculation of Eadm and 
P, the ADM mass is underestimated by about 0.1%, and 
the linear momentum is overestimated by about 0.3%. 
These errors cause the calculated radiation content to 
be too low. This is clearest for small values of P/M, 
for which the calculated maximum radiation content is 
negative. However, these are systematic errors, and do 
not prevent comparison between data sets generated and 
analyzed by the same procedure. 
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FIG. 8: The linear momentum P x of one CTSP black hole as 
a function of the value of (3q at the puncture. 




FIG. 9: The disagreement in linear momentum P x between 
the Laguna and full CTSP solutions for one black hole, as a 
function of the value of /3q at the puncture. 



shows the maximum radiation content of a number of 
single black-hole initial-data sets, and a comparison with 
similar cases for Bowen-York puncture black holes. The 
Bowen-York puncture data were generated by solving the 
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FIG. 10: Radiation content of single black-hole initial-data 
sets, with m — c = 1. Comparison of Laguna case with full 
CTSP data. 



FigurcHUIshows that single boosted CTSP and Bowen- 
York puncture black holes have the same maximum radi- 
ation content. They are indistinguishable from the point 
of view of the covariant quantities that we can calculate, 
the total ADM mass of the spacetime, the bare mass, and 
the linear momentum. 

There remains one free parameter in the CTSP ap- 
proach, the constant c in the lapse decomposition Q2ip. 
We investigate the effect of different choices of the con- 
stant c on CTSP boosted black-hole initial data in Figure 
rm We look at the maximum radiation content of initial- 
data sets for c = 0, 1.0, 2.0 and 5.0. The only limitation 
on c is that the lapse not become zero or negative at any 
point on the computational grid; if this is allowed to hap- 
pen, the numerical code will encounter division-by-zero 
errors when calculating the extrinsic curvature via 11911 . 
It is clear from Figure^Jthat the choice of c has no affect 
on the maximum radiation content of the initial data. 
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FIG. 11: Comparison of radiation content of single black hole 
initial-data sets, for various choices of the lapse parameter 
c > 0. The choice of c has a minimal effect on the physical 
properties of the initial data. 



maximum radiation content of single boosted Bowen- 
York black holes. This result suggests that single boosted 
CTSP black holes are nearly indistinguishable from single 
boosted Bowen-York black holes. It would be necessary 
to numerically evolve both types of initial data to verify 
this conjecture. These results are independent of the free 
slicing parameter, c. 

In the binary black-hole case, we do not expect CTSP 
and Bowen-York initial-data sets to be equivalent. Tichy 
and Briigmann 5] have shown that binary Bowen-York 
black holes do not satisfy = 0, and therefore cannot be 
a solution of the CTSP equations with that quasicquilib- 
rium choice. Similarly, there is no reason to expect that 
the physical properties of binary black-hole initial-data 
sets will be independent of the slicing parameters c,-. A 
detailed study of binary black-hole initial data will be 
made in a future publication. 



VII. CONCLUSIONS 

We have completed the incorporation of the conformal 
thin-sandwich decomposition of the initial-value equa- 
tions of general relativity into the puncture framework 
that was begun in [17|. We have developed a technique 
for constructing boosted black holes in the CTSP frame- 
work by specifying a value for the shift vector at the punc- 
tures, and have numerically constructed single black-hole 
initial-data sets using this prescription. 

We have studied the maximum radiation content of 
boosted CTSP black holes, and found that it equals the 
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